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We study the optical conductivity in the low-energy regime of gapped mono- and bilayer graphene. 
A scaling relation is found, in which the four parameters frequency, gap, Fermi energy and tem- 
perature appear only as combination of three independent parameters. The ratio of the optical 
conductivity of bilayer and mononlayer graphene is exactly 2. 

I. INTRODUCTION 

Monolayer and bilayer graphene are scmimetals with good conducting properties [1-3]. Moreover, 
in the presence of a microwave field the related optical conductivity is constant over a large regime of 
frequencies. This has been found in theoretical calculations [4-8] and was also observed experimentally 
[9, 10]. The reason for this behavior is the existence of at least two bands in both materials, where at 
Fermi energy Ep = (i.e. graphene without a gate potential) the lower band is occupied and the upper 
band is unoccupied. As consequence, the absorption of photons of energy Tiui from the microwave field 
creates electron- holes pairs due to the excitation of electrons from the lower band at energy — Too/ 2 to 
the unoccupied upper band at energy Tilo/2. This mechanism applies also to gated graphene which has 
a shifted Fermi energy Ep ^ 0. However, in this case photons can only be absorbed if Tilo/2 > Ep (for 
Ep > 0), since all the states in the upper band are occupied up to the energy Ep. Correspondingly, a 
photon can only be absorbed for —Tiw/2 < Ep if Ep < 0. This means that electron-hole creation by a 
microwave field is only possible if hui > \Ep\ 

It has been found in a number of recent experiments that the creation of a gap in the semimetallic 
band structure of monolayer graphene (MLG) is possible by absorption of hydrogen [11] or in bilayer 
graphene (BLG) by applying a double gate [12]. In both cases an electron-holes pair can also be created 
but this requires a photon energy larger than the band gap A (cf. Fig. 1). 

Once electron-holes pairs have been created they will contribute to a current in the material, where the 
latter is related to the strength of the external microwave field by the optical conductivity a^^(u). This 
quantity can be measured experimentally and characterizes the electronic properties of the material. In 
particular, it can be used to determine the band gap A, since it vanishes for Hlo < A. BLG, in contrast to 
MLG, has two low- and two high-energy bands. As a result, there are several gaps that lead to electron- 
hole pair creations on different energy scales with a more complex behavior of the optical conductivity 
[13-15]. 

In the following the optical conductivity shall be evaluated via the Kubo formalism for the low-energy 
bands in MLG and BLG at nonzero temperature T. This avoids electron-hole pair creation from higher 
energy bands and van Hove singularities. An important question in this context is the role of the low- 
energy quasiparticle spectrum on the optical conductivity. In order to focus on simple spectral properties, 
we consider only non-interacting electrons in a periodic tight-binding model. Thus disorder, electron- 
electron interaction and electron-phonon interaction are not taken into account. 

II. MODEL CALCULATION 



The low-energy quasiparticle states in MLG with a gap A = 2m are described by the massive two- 
dimensional Dirac equation 

H ( * t \ ( ~ E + m idx+d y 
M \yJ {id x -d y -E-m 

For simplicity, we have set the Fermi velocity vp = 1 because this parameter will not appear in the final 
results for the conductivity. A similar equation exists for the low-energy quasiparticle states of BLG 



5 



FIG. 1: Schematic picture of the creation of an electron-hole pair in gapped mono- or bilayer graphene by the 
absorption of a photon with energy Tkj. For this process the photon energy must be larger than the band gap 
A = 2m. 



[16, 17] 



With the plane-wave ansatz ^>i(x,y) = \J/^e lfca,x+ yV we obtain for MLG the following relations 



*d*,v) = 'Z,'''" V *<t(x,v), E 2 =m 2 + k 2 (3) 
E + m 



and for BLG 



* i (x,y)= (k * + ff ^(x,y), E 2 = m 2 + k i . (4) 
m + hi 

These solutions will be used as a basis for evaluating current matrix elements and the optical conductivity. 

Kubo formula: The optical conductivity can be calculated from the Kubo formula. This requires the 
evaluation of the current operator j M = ie[H, r^], where ry, (p = 1,2) is a component of the position 
of the quasiparticlc. The nonzero matrix elements of the current operator with respect to the energy 
eigenstates of Eqs. (3), (4) are either diagonal elements ie(Ek\[H,r,j]\Ek) or the off-diagonal elements 
ie(Ek | [H, r^]\ — Ek). It turns out that the diagonal elements do not appear in the real part of the optical 
conductivity [18]. Only the off-diagonal terms contribute because the optical conductivity requires a 
scattering process between two states whose energy difference is just the photon energy u>. A convenient 
representation of the Kubo formula then is [18] 

<r^Re M = ^-J | (E k | [H, r„] | - E k) \H(E k + „/2) = + ^ |* 

with the Dirac- Fermi distribution fp(x) = 1/(1 + e l3x ) at inverse temperature (3 = l/ksT. (Here and in 
the rest of this paper the spin and valley degeneracy, providing an extra factor 4, has not been written 
explicitly.) Integration over E k gives 



,2 , T ff>(E F + u/2) - U(E F - w/2) 



\(w/2\[H, ril }\-u>/2)\*dtpJ Jpy 1 ' Jpy ^ (5) 

with the Jacobian J from the radial k integration: 

v ' ; [ cj/2Vw 2 - 4to 2 BLG 
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current matrix elements: The commutator in the current operator is for Dirac fermions the Pauli 
matrix a^: [Hj/,^] = icr„ (MLG) and for BLG in Fourier representation 



[Hb,x] 



.9Hb 

' dk x 



2i(k x ai + k y <j 2 ) 



Then for the current matrix clement for MLG we obtain 



E 2 



which after angular integration yields 



\{E\a 1 \-E)\ 2 dp = n 



2m 2 



E 2 



m 2 +E 2 
E 2- 



For BLG (with k 2 = \Juj 2 /4 — m 2 , k x = k cos tp and k y = k sin ip) we have 



\{E\k x a 1 + k y a 2 \-E)f 



k 2 . 9 \ 2 -2 ^ ■ 2 

--^-+4— (cos i/j-sm (/jjcos ^ sin" 99 + — sin <p- 



and after the angular integration 



/ \{E\k x a x +k y <j 2 \ - E)\ 2 d(p = Tvt 
Jo 



,m 2 + E 2 
E 2 



2m 
~E 



COS ip Sill ip 



(6) 



(7) 



(8) 



(9) 



This is valid only for to 2 > 4m 2 = A 2 . As an example, these current matrix elements are plotted for 
E = 2 in Fig. 2 with and without gap. 

conductivity: Now we insert the results of Eq. (7) into the Kubo formula Eq. (5) and obtain for MLG 



a MLG 



8h 



A 2 ' 

1 + — 

u! 2 



&(uj 2 - A 2 )[fp(E F + w/2) - f (E F - «/2)] . 



(10) 



Inserting Eq. (9) into the Kubo formula gives exactly twice the conductivity of MLG: cr' BLG = 2a' MLG . 
Thus the conductivities of MLG and BLG agree up to a factor 2. The additive correction due to the gap 
parameter A 2 decays like lu~ 2 . As an example, the behavior of the conductivity versus uj/ksT is plotted 
in Fig. 3. 



III. DISCUSSION 



The conductivity in Eq. (10) obeys a scaling relation of the conductivity: 

a' xx (/3,A,E F ,uj) = a' xx (/3A,l3E F ,pLj) = a' xx (A/u,l3E F ,pu) , 

i.e. the conductivity depends only on three independent parameters. A similar scaling relation exists also 
for the conductivity in the presence of a scattering rate [5] . This implies that a reduction of temperature 
is equivalent to a simultaneous increase of frequency, gap and Fermi energy. 

The optical conductivity vanishes for photon energies less than the gap. It jumps to some finite value 
when the photon energy exceeds the gap energy, where the size of the jump depends on the Fermi energy 
(cf. Fig. 3). The maximal jump appears at E F = 0. This allows us to measure the gap by measuring 
the jump of the optical conductivity. Such a measurement can be performed as an optical transmission 
experiment [9, 19] because the optical conductivity is related to the transmittancc T through the relation 



1 + Stht^/c ■ 

The parameter c is the speed of light. Then graphene is completely transparent for photonic energies up 
to the gap energy. This could be used to filter light with frequencies higher than the one given by the 
gap energy. 
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FIG. 2: The current matrix element |{i5|[H, x]\ — E}\ 2 as a function of the angle tp between the wave vector k 
and the direction of the current operator. The left panel is the matrix element for E = 2 without gap (m = 0) 
and the right panel is the matrix element with gap (A = 2m — 2). Full (dashed) curves are for BLG (MLG). 
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FIG. 3: Optical conductivity of BLG as a function of ftui for /3m = 5 and /3Ef = 0, 4, 6, 10 with decreasing 
conductivity. These curves have to be multiplied by 1/2 to get the corresponding values for MLG. 



Our model assumption of taking into account only the low-energy bands of BLG restricts the photon 
energies Tilo to less than 0.8eV, which is the gap between the low-energy and the high-energy bands in 
BLG [13, 14]. This restriction also avoids stronger deviations from the Dirac theory of MLG and van 
Hove singularities. Taking into account high-energy bands does not change this picture qualitatively, 
since it would lead to additional jumps of the optical conductivities as soon as the photon energies exceed 
the gap energies. A van Hove singularity would appear as an additional peak. The effect of the gap is a 
global enhancement by to 2 , where for BLG the situation is more complex than for MLG. 

In conclusion, focusing on the low-energy dispersion of gapped monolayer and bilaycr graphenc, we 
have found simple expressions for the optical conductivities. They agree for all parameters up to a factor 
2. This is remarkable because the current-matrix elements of both graphene systems are rather different. 
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